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We show that, under conditions of electromagnetically induced transparency (EIT), a significant 
portion of the incident probe pulse can be transferred into Rayleigh and Raman scattering channels. 
The light scattered into the Rayleigh channel emerges from the sample with an EIT time delay. We 
show that a proper description of the probe light propagation in the sample should include, in the 
diffusion dynamics, a spin polariton generated by the two-photon EIT process. The results have 
important implications for studies of weak light localization, and for manipulation of single and few 
photon states in ultracold atomic gases. 
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I. INTRODUCTION 

External control of the electrodynamic response of 
atomic systems has been revolutionized by the merger of 
the ideas of coherent population trapping with 
the techniques and concepts of ultracold atom physics 
0- The potential afforded by the combination was first 
demonstrated in the remarkable experiments of Hau, et 
al. , in which a coherent light pulse was compressed in 
an ultracold gas of sodium atoms, the excitation having 
a very small group velocity ~ 17 m/s. Subsequent ex- 
tensive theoretical and experimental research has shown 
that a combined atomic-photonic quasiparticle excitation 
is created in the ultracold medium, and that the prop- 
erties of the polariton can be dynamically manipulated 
through the external qualities of the light fields used to 
prepare and probe the medium 0, 0; S ■ Since then, 
studies have shown a broad range of physics associated 
with coherent manipulation o f propag ation of electro- 
magnetic waves 11 la in in 111 111 Hi- For example, 
Chaneliere and co researchers have recently shown, in a 
series of beautiful experiments , that it is possible to 
generate single photon wave packets, and to map them 
into polaritons in an ultracold sample of ®^Rb atoms. 
The single photon wave packets could be regenerated af- 
ter a controllable delay by judicious application of a con- 
trol electromagnetic field. The physical processes associ- 
ated with manipulation and storage of individual photon 
wave packets, and entanglement of the quantum states 
of the packets with a propagating one, are critical ele- 
ments for quantum information protocols and quantum 
memory applications. 

Among the essential elements needed for practical ap- 
plications is quantitative understanding of coherence loss 
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mechanisms in each stage of a quantum information pro- 
tocol. In many studies of storage and retrieval of single 
and multiple quantum wave packets in ultracold gases, 
the initially created excitation undergoes decay as the 
length of the storage time is increased. In this paper, we 
focus quantitatively on one mechanism which can lead 
to such a potential loss of fidelity, and some of the sur- 
prising physics that results from these considerations. In 
particular, we consider the diffusely (multiply) scattered 
light that is necessarily generated as a result of optical 
excitation with temporally finite light pulses. In a gen- 
eral case, the associated finite spectral band width pro- 
duces diffusely scattered Rayleigh and Raman modes in 
the sample. These modes show quite different temporal 
behavior, in comparison with the coherent forward scat- 
tered probe light. In addition, because the Rayleigh scat- 
tered component maintains coherence with the incident 
probe beam, a new type of quasiparticle, a diffuse po- 
lariton, is created. Finally, we point out that the physics 
of light diffusion in a coherent medium consisting of a 
thermal (not ultracold) gas of Rb atoms has been inves- 
tigated by Matsko, et al. p^f. We point out that there 
have been a large number of studies of the fascinating 
physics of coherent population trapping; some represen- 
tative papers are cited here H 0, 0, Q| ■ A 
number of reviews have also appeared on the subjects of 
EIT and coherent population trapping, and contempo- 
rary applications [E II II Ellll ■ 

In the remainder of this paper, we first consider in 
more detail why, in the context of the previous discus- 
sion, the diffuse scattering channels are important. This 
is followed by a description and discussion of the complex 
susceptibility for the case of a typical lambda configura- 
tion attainable on hyperfine resonance transitions in an 
ultracold gas of *^Rb atoms |0| . The scattering channels 
are then considered, with particular attention paid to the 
Green's function for propagation in the inhomogencous 
and optically anisotropic medium. We then present our 
results describing the temporal behavior of the forward 
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scattered light and the diffusely scattered fluorescence. 
We finally show that the coherence of the multiply scat- 
tered Rayleigh mode can be detected through the appear- 
ance of the coherent backscattering effect for the weak 
probe beam. 

II. BASIC ASSUMPTIONS AND 
CALCULATION APPROACH 

A. Why the scattering channels are important 

As an example of a system where EIT resonance can 
be clearly observed, we consider the A-type configura- 
tion in the hyperfine manifold of the Di-line of ^^Rb, see 
Fig. ^ A strong coupling field with right-handed circu- 
lar polarization is applied between the F = 2 hyperfine 
sublevel of the ground state and F' = 1 hyperfine com- 
ponent of the excited level. The atoms equally populate 
the relevant Zeeman states of the lower hyperfine sublevel 
F = 1. In the figure, the detuning of the coupling (probe) 
laser from atomic resonance is defined as Ai (A2). Such 
an atomic configuration is typical for EIT observation 
and can be fashioned in ultracold systems using for ex- 
ample, magneto-optic or quasistatic dipole traps. Such 
configurations have also been widely studied in more tra- 
ditional heated atomic gas cells; see for examples, Ref. 

miiiiii. 

In the present paper, we assume that atoms are so 
deeply slowed by an atomic cooling process that we can 
neglect all effects associated with atomic motion. These 
are ideal conditions for EIT and the probe mode, applied 
with any polarization to the F — I F' = I hyperfine 
transition, would find respective A-type excitation chan- 
nels, thus generating, as a result of the coherent Raman 
process, hyperfine coherence in the ground levels. To set 
our terminology, we refer in brief to Rayleigh and elastic 
Raman scattering processes as Rayleigh scattering. The 
terminology Raman scattering refers only to hyperfine 
inelastic Raman scattering. All time scales are measured 
relative to the inverse lifetime of the ^^Rb resonance 
transition. 

To reliably observe the conversion of the probe pulse 
into a polariton-type quasiparticle state, which, as ex- 
pected, would forwardly propagate through the sample, 
some temporal and spectral requirements should be ful- 
filled. There are two important temporal parameters for 
the process: the pulse duration Tp and the delay time 
between the entrance and emergence of the probe light 
pulse. Apparently, for high-fidelity conversion of the full 
incident light pulse into the polariton pulse, it is desir- 
able that Td > Tp. For the simplest A-type configuration 
the delay time can be estimated as 

Td^ ^ ^ noX^L -^2" (2.1) 

where L is the sample length, w is a group velocity for the 
probe mode at the resonance, uq is the density of atoms. 
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FIG. 1; An example of an excitation scheme for observation 
of the EIT effect in the system of hyperfine and Zeeman sub- 
levels of the _Di-line of *^Rb. The coupling field is applied with 
right-handed circular polarization to the F — 2 F' = 1 
transition and the probe mode in the orthogonal left-handed 
polarization excites the atoms on the F = 1 F' = 1 tran- 
sition. The EIT effect appears for equal detunings of the 
coupling and probe modes from atomic resonances: Ai = A2 

X is the wavelength divided by 27r for the probe radiation, 
r is the natural decay rate for the upper state and flu 
is the Rabi frequency for the coupling mode. The time 
Tp can be estimated via the time-frequency uncertainty 
principle, i.e. as Tp ^ 1/Ap, where Ap is the spectral 
width of the pulse. The spectral width Ap is restricted 
by the condition that, at the relevant detuning from the 
EIT resonance, the optical thickness of the sample b{A2) 
at A2 = Ap would be small enough and the medium 
would be transparent. Then in order of magnitude the 
pulse duration Tp is limited by the inequality 



Tp > ^noX^L^ (2.2) 

R 

For an optically extended medium with \/ noX^L ^ 1 
the pulse duration can be made shorter than the delay 
time. But in reality such an optically dense sample is 
rather difficult to prepare in an experiment with ultracold 
atoms, where the parameter ngA^L is typically close to 
ten or even less in order of magnitude. 

The problem is more subtle if the EIT channel is ad- 
justed for transport of a portion of "non-classical" light. 
As a particular example, one can imagine a pulse of for- 
wardly propagating squeezed light. Such states of light 
can be created with an intra-cavity optical parametric 
light source and its spectral properties can be controlled 
with the quality factor of the cavity. For such a light 
source the outgoing radiation can be properly described 
by the model proposed by CoUett and Gardiner in Ref. 
[23 when the squeezed (Xi) and anti-squeezed (^2) 
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quadrature components are described by the correlation 
functions with different relaxation times, which we re- 
spectively denote as ti and T2. In principle, there exists 
an inequality between these correlation times ti > ti, 
which can be a strong constraint for high degrees of 
squeezing. For reliable transport of the squeezed light 
with preservation of its unique statistical properties it 
would be necessary that Tp > T2 > ti. The variation 
of all these temporal parameters should be limited both 
by the longer estimation (|2.1|) and the shorter estima- 
tion 1)2. 2|l . Apparently it would be more difficult to fulfill 
the EIT criteria for the squeezed light than for a pure 
coherent light source. 

However, the atomic sample can be probed with a light 
pulse with shorter duration than is given by (1221) and the 
EIT-mechanism will work even if a portion of the light 
pulse is transported via a non-forward scattering channel. 
For such an experimental situation the pulse duration 
Tp can be limited by the time scale r/f2|j, given by the 
inverse spectral width of the transparency window in the 
local susceptibility of the medium. Then the input pulse 
of the probe light should be transformed into a polariton- 
type pulse in the diffuse mode, which now will fill the 
sample via a coherent diffusion (non-forward scattering) 
process. 



where no(r) is the local density of atoms at a spatial 
point r. In a environment characteristic of ultracold and 
trapped atoms, the density distribution is typically in- 
homogeneous and the density matrix as well as the sus- 
ceptibility tensor are spatially dependent. The frequency 
detuning A2 is the offset of the probe mode uj2 from the 
resonance A2 = uj2 — (^f'f with F = 1, F' — 1. 

The second line in ()2.3|l reveals the contribution of 
the EIT-effect. Here Vnm' are the transition matrix el- 
ements for the coupling mode between those quantum 
states \n) = \F',n) and |m') = \F,m') which are subse- 
quently chained with an initial state |m) via the respec- 
tive A-type excitation channel. This is indicated in the 
sum by the dependence of the subscript indices on m: 
m! = m'ijn) and n = n{m). The frequency detuning Ai 
is the offset of the coupling mode cui from the resonance 
Ai = — lvf'f with F = 2,F' = 1. The pole in the de- 
nominator of Eq. H2.3|l is shifted due to the Autler-Townes 
effect and the self-energy correction is given by 

S„m'(A2) = A„m'(A2) — -r„m/(A2) 



B. The macroscopic susceptibility and scattering 
tensors 

The dynamics of the macroscopic polarization, induced 
by a probe radiation pulse, is driven by the dynamical 
susceptibility of the medium. Referring to the excita- 
tion scheme shown in Fig. ^ under conditions of EIT- 
resonance, the susceptibility of the medium for the probe 
mode entering the sample in a particular polarization is 
generated by two coherently interacting A-type channels. 
In the laboratory frame, with the z-axis directed along 
the coupling beam, the susceptibility tensor has a diago- 
nal form in the basis of circular polarizations and can be 
written as the following sum 



X^'(r,A2) = -<5,'' Yl 
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1 



h A2 + iT/2 



A2 + ir/2 Ai - A2 S„m'(A2) 



Pmm (^) 



(2.3) 



We use standard co/contravariant notation for the basis 
vectors of circular polarizations, see |21j], which can be 
expressed by Cartesian basis vectors as eo = e^, e±i = 

The first line in 1)2.3(1 has an isotropic form and de- 
scribes the local macroscopic susceptibility in the normal 
approach of linear electrodynamics. Here the squared 
transition dipole moments (de*)„m between the lower 
|m) = \F,m) and upper |n) = \F',n) Zeeman sub- 
levels are weighted with the population components of 
the atomic density matrix Pmm{^) = n'o(r)/{2F -1-1), 



The susceptibility tensor 1)2. 3|l describes an anisotropic 
and optically active medium despite the homogeneous 
population of the Zeeman suhlevels. This is a direct con- 
sequence of the EIT-effect. Anisotropy comes from the 
different A-type transitions activated by differently polar- 
ized probe modes. This effect has a similar physical na- 
ture as various optical anisotropy effects associated with 
the presence of the strong couphng field, see ^23, 24, 2bl. 

The scattering process in a medium is conveniently de- 
scribed by the scattering tensor formalism. This tensor 
is responsible for frequency- and polarization-dependent 
transformation of an incident electromagnetic plane wave 
as a result of its scattering on an isolated atom. Under 
conditions of EIT control for the mode incident on atom 
at frequency uj2, the scattering tensor is given by 

a(""™)(A2) ^ 4f ")(A2)|m")(m| 



y - 



m' (n},n 



1 idp^m" n{d'q)nra \ n\ i \ 

n A2 + ^r/2 '"^ 



X 1 



IK: 



A2-Kzr/2 Ai-A2 + S„„.,(A2) 



(2.5) 



which determines the amplitude of the outgoing wave for 
either the elastic or inelastic scattering channel accompa- 
nied by transition of the atom from the state |m) = \Fm) 
with F 1 to the state |m") = \Fm") with F=l 
(Rayleigh channel) or F = 2 (Inelastic Raman channel). 
The scattering tensor should obey the following impor- 
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tant identity 



47ra;2 



Im 



("™)(A2)(£*)Pe« 



(2.6) 



which is an optical theorem reflecting the unitary prop- 
erty of the scattering process. Here and (e'*)^ are re- 
spectively the contravariant circular components of the 
polarization vectors for the incident photon and of the 
complex conjugated polarization vector for the scattered 
photon. The sum is extended over all possible scattering 
channels and all possible output polarizations. After in- 
tegration over the full scattering angle the left hand side 
reproduces the total scattering cross section. 

Because of their similar physical nature, the expres- 
sions for the susceptibility tensor (|2.3|l . and for the scat- 
tering tensor H2.5|l are visualized in similar form. The 
important difference is that, in the expression 1)2. 5|l . an 
arbitrary reference frame can be assumed and any dipole- 
type transitions |m) —^ \n) can be initiated. At the same 
time, each of the upper states \n) is coupled by the strong 
field with a certain selected ground state \m') that is indi- 
cated in the sum as m' = m' [n). This prohibits the light 
scattering for the resonance mode when A2 = Ai but 
opens the scattering channel for any non-vanishing de- 
tuning. Moreover, under the EIT effect it becomes quite 
important which scattering channel (elastic Rayleigh or 
inelastic Raman) the probe photon occupies as it is scat- 
tered away from its original forward propagation. An 
elastically scattered photon still undergoes coherent cou- 
pling and should emerge from the sample with a rela- 
tively long delay. In contrast, in the Raman channel the 
photon freely propagates through the sample with nearly 
the (vacuum) speed of light because the medium is ideally 
transparent in this case. 



C. The Green's function 

Since the probe or scattered beam can propagate in 
any direction, its free-path transformation can be prop- 
erly described by the Green's function formalism. For 
the transmitted or elastically scattered light the Green's 
function can be expressed in the form of phase integrals 



perpendicular to the coupling beam. In this sense, ex- 
pression 1)2. 7|l gives a radiation zone asymptote for the 
retarded-type Green's function and the polarization in- 
dices qi , (72 are defined in the reference frame naturally 
linked with the ray direction, such that z \\ ri — r2. 
In this frame the Green's function has a non-diagonal 
form in a natural basis of local circular polarizations with 
qi^2 = ±1 (left-hand/right-hand). 
The phase integrals are given by 



0o(ri,r2) 



27rw2 



' r2 

ri,r2) = / x(r,A2)ds 2., 



I Xo(r,A2)ds 



and evaluated along a ray from point r2 to point ri. The 
integrand xo(r, A2) is expressed in terms of major-axes 
components of the susceptibility tensor as 



Xo(r,A2) = i±£25!^[x:l(r,A2)+xt}(r,A2) 



sin (3 Q 



x;](r,A2) 



(2.9) 



where (3 is the polar angle between the ray and the cou- 
pling beam directions. The function x(r, A2) is given by 



x(r,A2) 



Ex?(r,A2) 



Li=l 



1/2 



(2.10) 



where 

Xi(r,A2) = 

X2(r,A2) = 

X3(r,A2) = 



^cos7{i[x:J(r,A2)-fx:J(r,A2)] 



xg(r,A2) 



sin^/? 



2 sin7 -[xl}(r,A2)+xi}(r,A2)] 



cos/? 



xS(r,A2 



[x+J(r,A2)-x:J(r,A2)] (2.11) 



(ri,r2, 0^2) 



„i0o(ri ,r2)-|-jfc2 |ri -r2 



|ri -r2| 

X cos((/.(ri,r2))<55,''^ + isin((/)(ri,r2)) (n<?),,9^](2.7) 

Here k2 = 1^2 /c and n = fi{l^2) is the "unit" symbolic 
vector, the components of which are defined below. 

The parameters of the phase integrals can be expressed 
by the major components of the susceptibility tensor, 
which are defined by Ea. (|2.3|) and related to the probe 
beam either co-propagating in the forward direction or 



Here 7 is an azimuthal angle indicating the uncertainty in 
directions of a pair of orthogonal axes fixing the reference 
frame in the plane transverse to the ray. The "unit" 
symbolic vector n is defined by its components as 



ni = ni(A2) 



Xi(i-,A2) 
x(r,A2) 



(2.12) 



and is independent of r. It is described in the general 
case by the set of complex components (ni, 722, 713). 
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III. RESULTS 
A. The forward transmitted pulse 

The graphs in Figs. El and 01 show how the original 
pulse profile is modified after propagation through the 
sample in an EIT channel. The calculations were done 
for two different pulse shapes, these being Gaussian- type 
(Fig. ^ and rectangular-type (Fig. O , with pulse widths 
of Tp = 100 r~^. In our numerical simulations we consid- 
ered as an example ^''Rb atoms, when the coupling field 
and the probe light excitations were respectively applied 
at F = 2 ^ F' = 1 (empty) and F = 1 -> F' = 1 
(equally populated) hyperfine transitions, as shown in 
Fig. n The coupling field is always right-hand circularly 
polarized and the probe field can have varying polariza- 
tion, being either left- or right-hand circularly polarized. 
The calculations were based on attainable parameters of 
an ultracold atomic cloud confined in a magneto optic 
trap. For a Gaussian atom distribution in the trap, the 
weak-field optical depth, on resonance and through the 
center of the trap is given by bo = •\/27rnocroro, where 
no is the peak density and ctq is the normal (not mod- 
ified by the EIT effect) resonance cross section, tq is 
the radius of the cloud. The dependencies shown in 
Figs. 121 and (31 reproduce the pulse transformation af- 
ter its passing through such an atomic sample with the 
size ro = 0.25 cm and peak optical depth bo — 50. All the 
parameters are in accordance with the condition that the 
delay time is longer than the duration Tp of the orig- 
inal pulse. The inequality > Tp can be controlled by 
appropriate choice of the Rabi frequency for the coupling 
field flc = 2\Vnm'\, which we defined with respect to the 
m' = {F = 2, M = -1} ^ n ^ {F' = 1, M' = 0} hyper- 
fine Zeeman transition. For the dependencies shown in 
Figs. 121 and O the Rabi frequency is given by flc — OAT. 

For the case of a Gaussian profile the output pulse 
preserves the original shape but becomes more extended 
in time. The highest fidelity channel for reproduction 
of the original pulse is obtained when the coupling and 
probe fields are in orthogonal polarization states (blue 
curve in Fig|2Il. The loss of the pulse energy is only 10% 
for this case. But for the pulse transmitted through the 
sample in the right-handed polarization channel there is 
up to 40% energy loss. This is a clear manifestation of 
the scattering process, which cannot be ignored, as we 
argued by simple estimations in Section III Al 

The modification of the pulse profile manifests itself 
even more dramatically if the probe pulse has a rectan- 
gular shape. Due to dispersion the original pulse shape 
is completely transformed in the extended medium to a 
Gaussian form and the losses of energy are more signif- 
icant for this case. That effect is shown by the calcu- 
lations of the outgoing pulse profile presented in Fig. O 
for two different polarization channels. The calculations 
were done for the initial parameters similar to those used 
in the calculation of the dependencies for Fig. [21 For the 
outgoing pulses shown in Fig. Othe losses are up to 60% 
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FIG. 2: Intensity profiles for light pulses transmitted through 
a sample of *^Rb atoms via the EIT channels. The coupling 
field is applied on the empty F — 2 ^ F' — 1 transition 
with right-handed circular polarization and the probe field, 
exciting the atoms on the _F = 1 ^ _F' = 1 transition, is 
either left-handed polarized (blue, short-dash curve) or right- 
handed polarized (green, longer-dash curve). The black curve 
indicates the original Gaussian-type pulse with a duration of 

Tp = 100 r-^ 
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FIG. 3: Same as in Fig. |5]but for a rectangular input pulse 
profile. 

for the left-handed and up to 80% for the right-handed 
polarization channels. 



B. The scattered pulse 

The plots in Fig. 01 reproduce the time dependence 
of instantaneous intensity for the fractions of the light 
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FIG. 4: The intensity profiles for tlie portion of the light 
pulse scattered at 90° to the direction of the incident pulse. 
The curves represent the Rayleigh (Rl) channel (_F = 1 — > 
F' = 1 -» F = 1) and the Raman (Rm) channel (F = 1 ^ 
F' = 1 ^ F = 2), with the input polarization state, and 
polarization channel of the emerging light as indicated in the 
caption. In all cases, the observation channel corresponds to 
detection of light with left-hand helicity (/i-). 



pulse, originally incoming with Gaussian profile shown 
in Fig. 13 and scattered by the sample in a direction 
orthogonal to the incident probe beam. Note that we 
designate the polarization of the scattered light in terms 
of its helicity, which is a natural choice when considering 
the radiation propagating towards a detector. In Fig. 4, 
the curves indicate our results for the Rayleigh scatter- 
ing channel (F=1^F' = 1^F = 1) and the Raman 
channel [F — 1 ^ F' = 1 ^ F — 2). Let us recall here 
that for the sake of brevity we do not distinguish between 
elastic Raman and Rayleigh channels. The calculational 
parameters were chosen the same as for the dependencies 
of Figs. 121 and 13 The solid and dotted blue curves illus- 
trate the portion of the scattered light pulse incoming 
with left-handed circular polarization ((T_ ) and emerging 
the sample in the left-hand helicity polarization ) . In 
turn, the short-dashed and long-dashed green curves il- 
lustrate the portion of the scattered light pulse incoming 
with right-handed polarization (o'+) and emerging the 
sample again in the left-hand helicity polarization (ft.-). 
Thus, for the intensity profiles shown in Fig. ^ the out- 
going light is always considered in the left-hand helicity 
polarization channel. 

As follows from these dependencies, the scattered light 
emerges the sample with a rather long delay in units of 
r~^. This is an indication of the slow light phenomenon 
typically associated with the FIT effect, but for the scat- 
tered light. The portion of the input pulse, distorted 
in the scattering process, propagates through the sample 
faster and is more extended in time than the transmitted 
pulse. This can be explained by the spectral structure 
of the scattered pulse, which has a deficit of near reso- 



nance photons for which the EIT mechanism works ide- 
ally. For the scattered light pulse leaving the sample via 
a Rayleigh channel, the group velocity is greater than for 
the transmitted incident pulse. The photons created in 
the Raman channel leave the sample with the speed of 
light in vacuo because the medium is fully transparent 
for them. The spectral notch in the distribution of the 
scattered photons leads to the beating effect in the time 
behavior of the output light intensity, which is clearly 
seen for all the scattering channels presented in Fig. ^ 

To show the difference between the usual incoherent 
scattering and scattering in the environment supporting 
the EIT-effect, let us turn to the physical background 
of the process. If the medium is illuminated by a single 
mode probe field, which is precisely in the resonance with 
the coupling field such that wi = W2 + W21, where ^21 is 
the hyperfine splitting in the diagram of Fig. ^ the entire 
atom-field system is described by the following density 
matrix 



(3.1) 



j = l,2 



where dots denote the contribution of non-coupled state 
and the sum j — 1 , 2 is extended over two A schemes 
shown in figure ^ The wave functions contributing to 
the expansion (|3.1|l are given by 



(j) 



is) 



(j)2 



n, 



(j)2 



(j)2 



= |2)j 



X I Field) 



(3.2) 



where |l)j and |2)j are respectively the left and right 
atomic spin states coupled with the j-th A-type excita- 
tion channel and ilc and flp ^ arc the Rabi frequencies 
respectively for the coupling and for the probe fields. For 
the sake of simplicity we presume both the Rabi frequen- 
cies to be real parameters. It is a crucial point of the EIT 
effect that despite the fact that even though the atomic 
substate given in the square-brackets of (|3.2() is not an 
eigenstate of the atomic Hamiltonian, the full wavefunc- 
tion l^fj) is an eigenstate of the entire atom-field system 
if the field is in coherent state. Thus under conditions of 
perfect A-resonance the system is not able to scatter the 
light. 

However if the probe mode is not in exact two-photon 
resonance the situation becomes more complicated. In 
the single atom case and as a first approximation there 
will be coherent beats with exchange of a photon between 
coupling and probe modes of the field, which are initiated 
by low frequency phase oscillation of the atomic spin co- 
herence. In a macroscopic system these beats are collec- 
tivized and transformed to a polariton wave packet cre- 
ated in the sample by a pulse of the probe field enveloped 
by a set of near resonance spectral modes. For a short 
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distance the polariton wave propagates in the forward 
direction but it will be scattered at long distances if the 
medium is optically extended. However the probe field, 
even being scattered via Rayleigh channels, remains in a 
coherent state and can be also enveloped by the modes 
inside the EIT spectral window. Thus the polariton wave 
does not disappear as a result of Rayleigh scattering but 
transforms to a diffuse coherent mode. The quantum na- 
ture of such a diffuse polariton can be verified if the con- 
trol field is switched off and switched on again with delay. 
Then the propagation of the polariton wave in the for- 
ward direction and via Rayleigh scattering channel will 
be stopped and be regenerated again with a time delay. 
Thus the scattered light, as well as the transmitted light, 
can be stored in the medium and recovered on demand. 
That is a unique property of the light propagation under 
EIT conditions. 



C. The coherent backscattering process 

A clear indicator of the coherence in the Rayleigh scat- 
tering channel is the appearance of the coherent backscat- 
tering (CBS) process. This effect is a typical example 
of a "which path" interference fashioned by a disordered 
medium. For a portion of a light wave emerging the sam- 
ple in the backward direction there is constructive or, 
for some special cases, destructive interference [2^ be- 
tween the scattering amplitudes along any multiple scat- 
tering chain. The various manifestations of the CBS phe- 
nomenon in ultracold atomic systems [2^12^12^ . and un- 
der different physical conditions, has been the subject of 
many discussions during the last decade, see Kupriyanov, 
et al. for a review [s^l ■ Some future prospects in this field 
have been summarized by Havey and Kupriyanov "sT] . In 
the context of light scattering in an atomic environment 
supporting the EIT effect, the CBS process can be inter- 
preted as interference between two unknown reciprocal 
paths for the diffuse polaritonic wave propagating inside 
the sample via Rayleigh scattering as shown schemati- 
cally in Fig. 13 The specific nature of the EIT-resonance 
essentially modifies the time behavior of the interference 
component of the outgoing wave packet in comparison 
with how it is usually observed for a pulse-type excita- 
tion of an opaque atomic medium. 

For the scattering process developing in the standard 
conditions of elastic scattering and for a rather long ex- 
citation probe pulse the "ladder" (non-interference) and 
the " crossed" (interference) components have an approx- 
imately similar dependence on time. The important dif- 
ference in their time behavior is mainly observable in the 
transient stages of the excitation process, see Ref. [30|. 
Inside the pulse the relation between the instantaneous 
magnitudes of these components, usually expressed in 
terms of the enhancement factor for the CBS output in- 
tensity, has typically no time dependence. The enhance- 
ment factor approaches a stable and constant value inside 
the pulse. This value reproduces the enhancement factor 




FIG. 5; A schematic diagram explaining the coherent 
backscattering phenomenon. The polaritonic waves scattered 
in the medium follow pairs of reciprocal scattering paths indi- 
cated by the solid and dashed traces. The emergent electro- 
magnetic field in the backwards direction can show interfer- 
ence between the propagating modes following the two paths. 



for the steady state regime in conditions of a single mode 
excitation of the atomic sample with a resonant probe 
light. A crucial contrast for the CBS process observable 
in an environment of the EIT effect is that the resonant 
carrier mode of the light pulse cannot be scattered and 
the spectral profile for the pulse fraction, deflected by 
the medium from its original forward propagation, al- 
ways has a narrow gap near this mode. That results in 
a quite unusual time behavior for the interference over- 
lap of two reciprocal fragments of the backscattered wave 
packet. 

The dependencies of Figs. ElandlTlillustrate the tempo- 
ral behavior of the instantaneous intensity for the fraction 
of the probe pulse scattered in the backward direction 
via the Rayleigh channels (T_ — *■ /i_ and cr+ — > The 
outgoing intensity is mainly contributed by the " ladder" - 
type terms such that the interference contribution, shown 
in the magnified bottom parts of the graphs, gives actu- 
ally only a few percent of the total intensity. For conve- 
nience in these figures the gray curves indicate the time 
profiles for the incoming and transmitted pulses. Let 
us point out that these profiles slightly differ from those 
which are shown in Figs. 13 and 0] because the aperture 
of the probe beams are selected differently for the two 
calculations. In Fig. [21the cross section of the beam was 
less than the cross sectional area of the atomic cloud, 
which is a natural requirement for measurements in slow 
or stopped light experiments. Alternatively, for observ- 
ing the CBS effect from all atoms of the ensemble it is 
more natural to make the cross-section of the probe beam 
bigger than the cross sectional area of the atomic cloud. 
The important feature of the dependencies displayed in 
Figs. Eland[7|is that the "crossed" -type terms lead to an 
oscillating interference enhancement near the midpoint 
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FIG. 6: The temporal profiles for the fraction of the light 
pulse scattered in the backward direction for the a- — > h- 
Rayleigh scattering channel. The blue solid curve is the total 
intensity profile and the dotted blue curve is the interference 
contribution. The short- and long-dashed gray curves respec- 
tively indicate the input and the transmitted pulses. For these 
graphs the ordinate only qualitatively reproduces the relative 
magnitudes of the pulse intensities. 
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FIG. 7: Same as in Fig. |S| but for the a+ h- Rayleigh 
scattering channel. 



of the backscattered pulse profile. 

The time behavior of the interference effect can be 
properly described in terms of the enhancement factor 



Is{t)+lL{t)+Ic{t) 
Is{t)+lL{t) 



(3.3) 



where Is{t)+lL{t) is the contribution of single scattering 
and of the ladder terms in multiple scattering and Ic (t) 
is the contribution of the crossed terms. The dependence 
of j]{t) on time is illustrated in the graphs of Fig. |S1 
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FIG. 8: The enhancement factor for a+ —> /i_ (green dashed 
curve) and a_ h- (blue solid curve) scattering channels. 



It is an intriguing consequence of the calculation results 
that the behavior of the enhancement factor is generally 
quite complicated because of the nontrivial manifesta- 
tion of different spectral components of the input pulse 
in formation of the backscattered fraction of the outgoing 
pulse in different scattering orders. The most important 
seems an oscillating enhancement near the midpoint of 
the scattered pulse, which corresponds with the polari- 
tonic wave packet overlap inside the medium. 



IV. CONCLUSIONS 

In summary, we have made a theoretical and calcula- 
tional investigation of the Rayleigh and Raman scattered 
light modes in a medium configured to demonstrate elec- 
tromagnetically induced transparency in the coherently 
forward scattered light. In particular, for an anisotropic 
and inhomogeneous sample of ultracold atoms we have 
examined the time evolution of the forward scattered 
light, the diffusely (multiply) scattered light in a char- 
acteristic right-angle fluorescence geometry, and in the 
backscattering configuration used to study weak local- 
ization of light in ultracold atomic gases. The inten- 
sity of the forward scattered light is modified by both 
Rayleigh and Raman scattering of Fourier components 
of the incident pulse which lie outside the EIT trans- 
mission window. In particular, the intensity of a light 
pulse regenerated in a typical 'stopped-light' experiment 
is significantly modified by the loss of these components, 
which in part maintain their coherence with respect to 
the incident light. In addition, these components pro- 
duce complex time dependence in both the fluorescence 
and in the coherent backscattering geometries. The in- 
terferences responsible for the coherent backscattering 
enhancement are a clear indicator of the coherence of 
the scattering processes. We emphasize the fundamental 
conclusion that light scattered into the Rayleigh modes 
maintains its coherence with respect to the original in- 
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cidcnt pulse. This is manifested by interfcrometric en- 
hancement of the backwards scattered intensity from the 
sample. This in turn implies that, along with the polari- 
ton associated with the forward scattered light, a diffuse 
and coherently related polariton coexists within the vol- 
ume of the ultracold atomic sample. We close by noting 
that a deeper understanding of the mesoscopic proper- 
ties of the diffusive excitation is clearly required, and 
this report represents a first step in this direction. Fur- 
ther theoretical and experimental research into this area 
is underway. 
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